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Cl  (Optical )  Studies  in  the  Free  Atmosphere 
Based  on  Rawinsonde  Data 


1.  INTRODUCTION 

Atmospheric  turbulence  has  Long  been  known  by  astronomers  to  affect  light 

radiation  passage  through  the  atmosphere.  For  optical  systems,  the  single  most 

2 

important  parameter  describing  the  turbulent  atmosphere  is  Cn  (the  index  of  refrac¬ 
tion  structure  parameter).  This  parameter  is  also  related  to  the  temperature 

2  2 
structure  parameter  (C„),  which  is  a  variable  that  is  easier  to  measure  than  C  . 

2  T  2  n 
Cn  can  be  obtained  from  measurements  of  C^.  using  the  simplifying  assumption 

that  the  index  of  refraction  fluctuation  depends  only  on  the  fluctuations  in  tempera¬ 
ture.  In  this  case  the  only  other  variables  required  are  the  temperature  and 
pressure.  The  temperature  structure  function  can  be  obtained  by  using  thermo¬ 
sondes  that  require  two  similar  probes  spaced  a  distance  r  apart  for  measuring 
temperature.  Thermosonde  measurements  can  be  obtained  "in  situ"  in  the  free 

atmosphere  using  aircraft  and  balloons.  Remote  sensing  techniques  for  indirect 
2 

measurements  of  C  include  stellar  scintillometers  and  incoherent  radars.  The 
n 

scintillometers  measure  the  fluctuations  in  starlight  caused  by  the  atmosphere  for 
different  spatial  wavelengths  and  use  weighting  functions  for  altitude  resolution. 
Radars  detect  returns  from  fluctuations  in  the  refractive  index  of  the  atmosphere 
caused  by  turbulence.  The  reader  is  referred  to  a  section  in  the  Infrared  Handbook 

(Received  for  publication  19  April  1984) 


entitled  "Propagation  Through  Atmospheric  Turbulence  for  more  detailed  back- 

2 

ground  information  on  C  and  other  optical  parameters.  This  section  presents 

n  2 

excellent  practical  applications  for  obtaining  estimates  of  C  and  also  contains  a 

n2 

large  number  of  references  on  the  subject.  Also,  Good  et  al  contains  compari- 

2 

sons  of  recent  scintillometer  and  radar  measurements  of  C  . 

n  3 

In  this  paper  several  models  are  reviewed  and  the  Hufnagel  model  is  compared 
4  5 

with  the  Van  Zandt  model.  *  Both  models  use  rawinsonde  thermodynamic  data 

2 

to  derive  estimates  of  C  .  The  rawinsonde  data  base  has  been  previously  used 

6  n  7 

by  Murphy  et  al,  and  Murphy  and  Scharr  to  model  Richardson-number  ( Ri> - 
inferred  estimates  of  the  likelihood  of  occurrence  of  turbulence  in  the  lower  atmos¬ 
phere.  The  same  derived  quantities  used  in  that  study,  namely  the  winds,  N  (the 
square  of  the  Brunt -Vaisalla  frequency,  S  (the  vertical  shear  of  horizontal  winds 

squared  term),  and  the  ratio  N/S  (the  Richardson  number)  are  required  and  used 

2 

as  inputs  to  the  Hufnagel  and  Van  Zandt  models  to  obtain  C  estimates.  Also, 

2  n  2 
drop-off  rates  of  Cn  and  coherence  lengths  (rQ)  are  determined  and  the  Cn  drop-off 

rates  are  compared  to  those  obtained  from  radar  measurements. 

The  twice  daily  rawinsonde  balloon  measurements  of  wind  and  temperature  at 

each  of  many  locations  (81  locations  in  the  continental  United  States)  are  fit  using 

a  function  that  closely  approximates  a  linear  fit  from  point  to  point  of  the  original 

data.  The  temperature,  the  temperature  gradient,  and  the  wind  component  gradients 

are  determined  from  the  fitting  function  equation  by  computer  at  each  kilometer 

level  from  1  km  to  25  km.  Daily  values  of  the  stability  term,  shear  term,  and 

Richardson  number  are  accumulated  in  bins  at  each  km  level  to  provide  seasonal 

statistics  <~  180  values  /season).  There  are  two  reasons  for  not  using  data  below 

1  km.  The  first  is  that  the  stations  used  in  this  study  had  varying  site  elevations 

that  precluded  obtaining  a  good  statistical  sample  for  heights  below  1  km.  Also, 

data  were  not  used  below  1  km  since  we  were  inclined  to  limit  this  study  to  that 

region  regarded  as  above  the  "boundary  layer"  which  is  generally  defined  to  be  above 

approximately  1  km.  The  upper  level  of  25  km  was  chosen  since  this  criterion 

provided  a  near  maximum  of  100  samples  per  season  for  all  stations  at  this  level. 

Balloon  bursts  above  this  level  greatly  reduce  the  sample  size. 

The  results  presented  here  represent  only  a  first  step  in  developing  a  climatology 
2  2 

of  Cn  for  optical  effects.  Cn  estimates  are  obtained  only  above  the  first  kilometer 

and  for  dry  conditions.  For  optical  effects  then,  we  are  only  interested  in  the  optical 
2 

refraotivity,  C  (opt),  which  is  the  radio  refractivity  less  the  humidity  term.  Also, 
n  2 

results  from  models  should  be  compared  with  actual  measurements  to  determine 
differences  which  have  to  be  accounted  for  by  improvements  in  the  models. 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  37.  ) 


2.  A  SURVEY  OF  SEVERAL  MODELS  FOR  C? 

2 

In  most  optical  applications,  one  must  integrate  the  product  of  C  and  some 

1  n 

function  of  distance  along  the  optical  path.  For  other  than  horizontal  propagation, 

2 

we  need  to  be  able  to  model  the  height  dependence  of  Cn  since  it  is  highly  unlikely 

that  values  will  be  directly  available.  Furthermore,  to  construct  climatologies 
2 

of  Cn  from  rawinsonde  profiles  of  wind  velocities,  temperature,  and  pressure,  a 

theoretical  model  is  needed.  In  this  section  we  discuss  several  such  models. 

2  1 

An  extremely  simple  model  for  Cn  is: 


below  20  km  above  sea  level 


above  20  km  above  sea  level 


<  1) 


where  z  (in  km)  is  altitude  above  local  ground.  The  advantage  of  this  model  is  its 

simplicity,  but  it  cannot  be  expected  to  yield  anything  other  than  qualitative  results 

when  used  to  calculate  the  effects  of  turbulence  on  optical  beams. 

1  3 

Hufnagel  '  used 


C2  2.7  X  ia'8 9 lo|3  u2(z/  10)10  e'z  +  e'z/ 1,51  . 
n 


In  this  empirical  model,  z  is  the  elevation  above  sea  level  in  km,  and  u  is  the  rms 
wind  speed  (in  m/sec)  in  the  range  from  z  =  5  km  to  z  =  20  km,  that  is, 

9  20  2 

u“  -  1/ 15  /  u  tz )  dz  (3) 

5 


where  u(z)  is  the  wind  speed  at  altitude  z.  This  model  was  constructed  to  best  fit 
2 

measurements  of  C  from  z  »  24  km  down  to  the  first  strong  inversion  layer  and 

n  2 
so  this  is  its  expected  range  of  validity.  Use  of  the  appropriate  value  of  u  ob¬ 
tained  from  rawinsonde  data,  brings  this  model  to  within  a  factor  of  2  of  the  median 
2  8  9  10 

values  of  Cn  obtained  by  other  methods.  ’  ’  Although  this  model  cannot  be 

8.  Balsley,  B.  B.  ,  and  Peterson,  V.  L.  (1981)  Doppler  radar  measurements  of 

clear  air  turbulence  at  1290  MHz,  J.  Appl.  Meteorol.  ,  20:266-274. 

9.  Barletti,  H.  .  Cappatelli,  G. ,  Paterno,  L. .  Righini,  A.,  and  Speroni,  M. 

(197  6)  Mean  vertical  profile  of  atmospheric  turbulence  relevant  for 
astronomical  seeing,  J.  Opt.  Soc.  Am.,  66:1380-1383. 

—  2 

10.  Nastrom,  G.  D. ,  Gage,  K.S.  ,  and  Balsley,  B.  B.  (1982)  Variability  of  C  at 
Poker  Flat,  Alaska,  from  mesosphere,  stratosphere  and  troposphere 
(MST)  Doppler  radar  observations.  Opt.  Eng.  ,  21:347  -351. 


expected  to  agree  with  a  single  set  of  measurements  of  C  ,  it  can  be  modified  to 

3  n  2 

yield  typical  fine  structure  patterns.  We  note  that  for  z  and  u  large  enough,  this 

2  2 
model  predicts  that  CR  is  directly  proportional  to  u  . 

Two  recent  theoretical  models  have  been  developed  by  Van  Zandt  et  al*’  ®  to 
2 

allow  calculations  of  C  from  routine  rawinsonde  data.  These  models  are  based 
n 

on  the  empirical  fact  that  turbulence  occurs  in  thin  layers  and  on  the  theoretical 
work  of  Tatarskii**  who  has  shown  that  in  a  homogeneous  isotropic  turbulent  layer 

Jl  2  -  _  4/3  .,2  ... 

C  =  a  a’  L  M  ,  (4) 

n 

2 

where  a  is  a  constant  approximately  equal  to  2.8,  a’  is  the  ratio  of  eddy  diffu- 
sivities,  which  is  assumed  to  be  unity,  L  is  the  outer  scale  length  of  the  turbulence 
spectrum,  and  M  is  the  vertical  gradient  of  the  refractive  index  (neglecting  humidity). 
M  is  given  by: 

M  =  -  77.  6  X  10"6  P/T  In  5  ,  (5) 

where  P  is  the  atmospheric  pressure  in  millibars,  T  is  the  absolute  temperature 
and  9  is  the  potential  temperature. 

2 

It  is  important  to  emphasize  that  this  relationship  for  holds  only  in  the  pres¬ 
ence  of  shear  turbulence.  This  is  defined  to  occur  when  the  Richardson  number, 

Ri,  is  equal  to  or  less  than  one -quarter,  where 

Ri  =  N/S  (6) 

andN=g^-ln0,  g  =  9,8  m/sec^,  and  S  =  (9V/ 9z)2.  N  is  the  static  stability  term 
and  S  is  the  shear  parameter.  (Note:  this  notation  has  by  no  means  been  standard¬ 
ized  in  the  literature  and  care  must  be  used  in  comparing  formulas  from  different 
sources.  We  use  N  for  the  square  of  the  Brunt -Vaisalla  frequency  and  S  for  the 
square  of  the  vertical  velocity  gradient.  )  On  substitution  we  obtain 


1.76  X  10 


10  (P/T)2  N2  L4/3,  if  Ri  =  N/S  %  1/4  . 


11.  Tatarskii,  V.  I.  (197  1)  The  Effects  of  the  Turbulent  Atmosphere  on  Wave 
Propagation,  U.  S.  Dept,  of  Commerce,  National  Technical  Informatior 
Service,  Springfield,  VA,  pp.  74-7  6. 


2 

This  result  shows  the  basic  problem  in  using  rawinsonde  data  to  calculate  Cn< 

What  is  obtained  from  this  data  base  are  values  of  N  and  S  averaged  over  hundreds 

of  meters,  whereas  the  use  of  the  above  equation  requires  the  use  of  local  values 

of  N  and  S  only  for  N/S  =  1/4,  In  the  first  model  of  Van  Zandt  et  al,  *  it  is  assumed 
2  2 

that  C  ,  the  average  value  of  C  across  an  observed  layer  of  atmosphere,  can  be 
n  ^2  ^2 

calculated  from  Cn  =  C*  F  where  C*  is  calculated  using  averaged  values  of  P,  T, 

and  N  for  the  given  layer  and  F  is  the  fraction  of  the  layer  that  is  turbulent,  L  is 

treated  as  an  adjustable  parameter.  F  is  calculated  by  assuming  that  the  probability 

distribution  for  the  velocity  shear  is  normal,  with  the  mean  value  obtained  from 

rawinsonde  data,  and  has  a  standard  deviation  of  0.  010  sec  *  in  the  troposphere 

and  0.  015  sec  *  in  the  stratosphere.  F  can  then  be  calculated  by  integrating  this 

distribution  over  those  values  for  which  Ri  t  1/4.  With  a  choice  of  L  =  10  m, 

2  2 

good  agreement  between  and  an  average  of  from  radar  data  is  obtained. 

The  above  model  has  been  improved  recently  by  Van  Zandt  et  al^  by  including 
fine  structure  in  the  static  stability,  N,  in  addition  to  the  velocity  shear  term,  S. 
Furthermore,  the  parameter  L  is  treated  as  a  variable  layer  thickness.  It  is 
assumed  that 

OO 

C2  =  /  dL  /  /  dS  dN  pl(L,  S,  N)  C2  (8) 

n  n  n 

°  R.  S  1/4 

where  p^L,  S,  N)  dL  dS  dN  is  the  probability  that  a  point  in  a  turbulent  layer  has 
values  of  L,  S  and  N  lying  between  L  and  L  +  dL,  S  and  S  f  dS  and  N  and  N  +  dN.  It 
is  assumed  within  the  limits  of  integration  that  p4(L,  S,  N)  =  0  if  Ri  >  1/4.  Since 
pl(L,  S,  N)  is  not  known  precisely,  certain  assumptions  are  required.  It  is  assumed 
that 

pl(L,  S.  N)  =  pt(L)pt(S)pt(N)  .  (9) 

that  is,  the  marginal  distributions  are  independent  and  further  that 

pl(X)  =  p°(X)  ,  (10) 

where  the  superscript  of  zero  indicates  non-turbulent  conditions  and  X  represents 
any  of  the  three  variables.  This  means  we  are  assuming  that  the  probability 
distributions  for  the  individual  variables  are  unchanged  as  we  change  from  turbulent 
to  non-turbulent  conditions.  p^(L),  however,  must  be  related  to  the  probability 
that  a  point  is  in  a  layer  of  thickness  L  whether  that  layer  is  turbulent  or  not.  The 
individual  distributions  assumed  are  as  follows: 


(a)  p  (S)  =  exp(-(S*  S)/2a“]  Ifll(SXS)  '  /as)/2os 


This  is  the  Rice-Nakagami  distribution.  "5  can  be  calculated  from  the  rawinsonde 
data  and  is  the  standard  deviation  of  either  component  of  8V/3z.  For  computa¬ 
tional  purposes,  the  Rice-Nakagami  distribution  is  rewritten  in  the  form 

P°(S)  =  exp  -(— S----)  exp  (-x)  1  (x)/2a^  (12) 


where 


x  -  (ssr)0-5/a^ 


Based  on  Essenwanger's  empirical  work  on  standard  deviations  for  wind  shears, 
it  is  assumed  that 

-1  -1 

ag(S  )  as<S  )L  s  (14) 

-2-1  -2-1 

with  as  chosen  as  2  X  10  sec  and  3  X  10  sec  for  the  troposphere  and 
stratosphere  respectively  and  a _  chosen  as  0.  3.  We  note  that  with  the  choice  of 
L  =  10  m,  this  is  consistent  with  the  197  8  model4  values  of  Og.  1Q  is  the  modified 
Bessel  function  of  the  first  kind. 

(b)  p°(N)  =  exp[-(N  -  W)2/20^]/<2tt)1/2  <?n  (15) 

1/2 

where  N*  is  obtained  from  rawinsonde  data  and  it  is  assumed  that  o  ^  =  (N)  1  <Tg 

(c)  p°(L)  =  exp(-L/D  [<L/D  (1/D]  .  (16) 

The  justification  for  this  distribution  is  based  on  the  likelihood  of  an  exponential 
distribution  of  L  with  the  factor  L/H  added  to  provide  zero  probability  at  zero  length. 
U  is  chosen  to  best  fit  the  radar  data  and  again  turns  out  to  be  10  m. 

These  assumptions  yield  better  agreement  with  radar  results. 


12.  Essenwanger,  O.  (1963)  On  the  derivation  of  frequency  distributions  of  vector 
wind  shear  values  for  small  shear  intervals,  Geophysics  Pura  e  Applicata, 
^6^:2  16-224. 


In  either  of  Van  Zandt's  models  we  can  see  the  somewhat  complicated  de- 

pendence  of  Cx  on  the  Richardson  number.  As  Ri  decreases,  the  fraction  of  the 

observed  layer  that  is  mechanically  turbulent  increases,  which  of  itself  would 

cause  C2  to  increase.  On  the  other  hand,  if  Ri  decreases  due  to  a  decrease  in  N 
n  '  2 

(the  static  stability)  this  would,  of  itself,  cause  C  to  decrease  since,  from  Tatarskii's 

2  ^2 

basic  result,  we  can  see  that  Cn  is  directly  proportional  to  N  . 

The  subroutine  shown  in  Appendix  A  was  used  to  implement  Van  Zandt's  1981 

13 

model,  with  two  changes,  suggested  by  Van  Zandt,  from  its  original  formulation. 

These  are  discussed  below.  The  inputs  to  this  subroutine  are  77  >  0  and  i?  obtained 

— T  2 

from  our  fitting  routine  (Section  1)  and  the  output  is  C  /(P/T)  ,  which  is  then  con- 

2  n 
verted  to  C  .  If  IT  <  0,  we  assume  that  the  entire  region  is  turbulent  and  calculate 

C2  from  Van  Zandt's  1978  model  with  F  =  1. 

n  -2-1 

The  form  for  a„  used  has  been  changed  from  2  X  10  sec  (troposphere)  and 

-2  - 1  ^  1/4 

3  X  10  sec'  (stratosphere)  to  ac  =  0.2(77)  .  This  one-fourth  power  dependence 

b  14 

on  N  is  supported  by  some  recent  theoretical  work  of  Weinstock.  The  advantage 
of  this  form  is  that  it  provides  a  continuous  variation  of  Og  and  does  not  require 
the  location  of  the  tropopause. 

The  distribution  function  p**(L)  has  been  changed  from  (L fl?)  e  to  a  flat 

distribution.  We  assume  p**(L)  is  0.  01  over  the  range  of  L  values  used  in  the  inte¬ 
gration.  Detailed  analysis  of  the  contribution  of  different  layer  thicknesses  to  the 
2" 

total  C  value  reveals  that  thick  layers  beyond  100  m  contribute  small  amounts. 

n  5 

This  definition  of  L  is  used  in  the  Van  Zandt  model  for  ease  of  computation.  The 

Tatarskii**  definition  of  L  is  the  outer  scale  length  of  isotropic  turbulence.  The 

definition  of  shear  turbulence  through  the  Richardson  number  criterion  leads  to  a 

turbulence  layer  thickness.  H.  used  to  define  the  vertical  gradients  of  potential 

temperature  and  horizontal  winds.  The  relationship  between  outer  scale  length,  L, 

15 

and  turbulent  layer  thickness,  H,  is  not  well  defined.  Experiments  by  Good  et  al 

have  indicated  that  above  the  boundary  layer,  the  outer  length  scale  is  about  1/10 

5 

the  turbulent  layer  thickness.  If  this  is  incorporated  into  the  Van  Zandt  Model, 
then  substituting  L  =  H/10  into  Eq.  (7). 

C2  ^  8.  2  X  10'*2  (P/T)2  N2  H4/3  (17) 

n 

2 

with  N  and  Ri  evaluated  on  the  basis  of  turbulent  layer  thickness  H.  With  this 
modification,  it  will  be  necessary  to  increase  the  probability  of  turbulent  layer 


13.  Van  Zandt,  T.  E.  ( 1983 )  private  communication. 

14.  Weinstock,  J.  ( 1980)  Vertical  wind  shears,  turbulence,  and  non -turbulence 

in  the  troposphere  and  stratosphere.  Geophys.  Res.  Letts. 

15.  Good,  R.  E.  ,  Quesada,  A.  F 


7:749  -752. 

W* 


Brown,  .).  H.  ,  and  Dewan,  E.  M.  (1984) 
Probability  distribution  of  turbulence  layer  thickness  in  the  stratosphere, 
■1.  Geophvs.  Res,  ,  to  be  published. 


occurrence  to  match  observed  optical  turbulence  values.  The  effect  is  to  change 

2 

constants  but  does  not  change  the  underlying  relationship  of  to  meteorological 
wind  and  temperature  data. 

3.  STATISTICAL  RESULTS  FOR  C*  FROM  VAN  ZANDTS  MODEL 

■ 

2 

Our  examination  of  the  characteristics  of  the  C  distribution  included  four 

n 

stations:  Barrow,  Alaska  (7  1°  N),  International  Falls,  Minnesota  (48°  N), 

Pt.  Mugu,  California  (34°  N),  and  Brownsville,  Texas  (25°  N)  for  each  season: 
winter,  spring,  summer,  fall.  Each  season  was  defined  using  solstice  and  equinox 
dates: 

Winter  -  December  21  through  March  20, 

Spring  -  March  21  through  June  20, 

Summer  -  June  21  through  September  20, 

Fall  -  September  21  through  December  20. 

9 

C  median  seasonal  values  were  computed  at  altitudes  ranging  from  1  to  25  km 
using  the  N  and  S  values  determined  from  rawinsonde  daily  measurements. 

A  series  of  plots  were  drawn  examining  seasonal  and  latitudinal  variations  of 

2 

C  as  a  function  of  altitude.  Different  combinations  of  stations  and  seasons  were 

n  2 
overlayed  to  examine  the  structure  of  the  median  C  values  over  1-25  km.  Fig¬ 
ure  1  contains  the  two  extremes  in  latitude  for  the  winter,  Brownsville,  Tex.  , 

(25°  54'  N)  and  Barrow,  Alaska  (7  1°  18'  N).  There  are  significant  differences  in 

2 

the  distribution  of  median  values  between  10  and  20  km.  Although  not  as  pro¬ 
nounced,  the  same  structure  was  found  between  Pt.  Mugu,  Calif.  (34°  6'  N)  and 
International  Falls,  Minn.  (45°  34'  N),  which  is  shown  in  Figure  2.  All  four 

stations  are  plotted  in  Figure  3.  As  latitude  increases,  there  is  a  distinct  decrease 

2  2 

in  C  between  10  and  20  km.  Examination  of  summer  median  C  latitudinal  varia- 
n  n 

tions  revealed  only  slight  differences  over  the  1  -  25  km  region  (Figure  4).  There 
was  a  shift,  however,  in  each  distribution  in  the  approximate  height  range  between 
12  and  15  km.  This  may  be  due  to  the  jet  stream. 


2 

A  comparison  was  made  of  winter  and  summer  median  C  values  for  each  of 
the  four  stations.  The  winter  and  summer  plots  for  Barrow,  Alaska  are  presented 
in  Figure  5.  The  winter  curve  was  fairly  linear  from  1  to  25  km,  while  the  summer 
curve  had  an  inflection  point  around  10  km.  International  Falls,  Minn.  (Figure  6) 
and  Pt.  Mugu,  Calif.  (Figure  7)  demonstrated  a  similar  trend  with  inflection  points 
at  10  and  13  km  respectively.  The  graph  of  Brownsville,  Tex.  (Figure  8)  did  not 
reveal  the  same  structure  as  the  others. 

2 

We  were  interested  in  seasonal  and  latitudinal  variations  of  C  for  each  of 

n 

the  stations  mentioned  above.  A  non -parametric  test  was  selected  to  measure 

the  distributional  differences  between  the  four  groups.  This  is  called  a  "median 
1 6 

test"  because  the  data  from  each  group  are  arranged  from  the  lowest  to  highest 

values  and  each  value  is  then  compared  with  the  median  of  the  combined  groups. 

If  the  four  groups  are  in  fact  from  the  same  distribution,  the  number  of  values 

above  and  below  the  overall  median  should  be  about  the  same.  The  chi-square 

statistic  measures  the  degree  to  which  each  group  deviates  from  the  expected 

combined  distribution.  A  contingency  table  is  then  arranged  and  analyzed  using  a 
2 

chi-square  (x  )  statistic  having  three  degrees  of  freedom  (number  of  groups  minus  1). 

o 

The  hypothesis  of  identical  distributions  is  rejected  if  the  observed  x  is  significantly 
large. 

There  was  no  significance  found  in  accepting  the  hypothesis  of  identical  dis¬ 
tributions  in  examining  latitudinal  variations  from  station  to  station  for  each 
season.  Each  of  the  altitude  levels  from  i  to  25  km  were  examined  separately 

(Table  1).  A  significantly  large  x2  statistic  at  the  5  percent  level  of  significance 
2 

was  greater  than  \  critical  =  7.  18.  As  is  apparent  from  Table  1,  random  altitude 
levels  demonstrate  similar  distributional  characteristics;  there  is  no  apparent 
physical  reason,  however,  why  these  exist.  Therefore  we  can  conclude  that  there 
are  distinct  differences  in  the  distributions  of  the  four  stations,  having  differing 
latitudes,  over  the  four  seasons  defined  above. 

In  the  examinations  of  seasonal  variations,  each  station  was  tested  for  differ¬ 
ences  season  to  season  at  each  kilometer  level  from  1  to  25  km  using  the  test 
previously  discussed.  A  similar  phenomenon  was  found  as  with  the  latitudinal 
variation  analysis  (Table  2).  There  were  significant  differences  season  to  season 
throughout  each  altitude  bin  from  1  to  25  km. 

16.  Dixon,  W.  J.  ,  and  Massev,  F.  J.  (1969)  An  Introduction  to  Statistical  Analysis 
McGraw-Hill,  p.  351.  ~ 
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Table  1.  C  Median  Test  Results  on  Station-to-Station  Differences 
n 


Altitude 

Winter 

Spring 

Summer 

Fall 

1 

16.77 

13.68 

22.  77 

$ 

6.93 

13.20 

❖ 

* 

2 

3.  96 

13.  96 

4.75 

3 

8.80 

9.  27 

25.  28 

0.  94 

4 

26.75 

12.  11 

52.  09 

2.  73* 

5 

25.25 

49.  59 

$ 

6.  29 

23.47 

6 

12.69 

39.42 

11.  78 

13.  13 

7 

6.  59* 

35.  20 

11.  08 

9.41 

8 

9.  21 

39.  16 

7.  64* 

18.  59 

9 

9.  50 

28.  38 

9.63 

15.41 

10 

18.  83 

28.  20 

22.60 

* 

7.41 

11 

82.63 

43.  92 

94.  29 

75.33 

12 

140.  19 

90.54 

26.70 

136.30 

13 

213. 70 

160.  14 

70.  33 

168.  21 

14 

229.  28 

230.  14 

171.  19 

216.50 

15 

242.40 

363.  35 

209.  91 

243.52 

16 

270.  83 

404.65 

260.  06 

345.  92 

17 

261.  89 

388.  56 

221.  85 

312.  08 

18 

188.40 

287.  23 

193.  36 

169.  55 

19 

154. 28 

262.  32 

163.49 

132.  57 

20 

102.  97 

174.  87 

140.62 

109.  73 

21 

44.  71 

117. 13 

105.  86 

46.64 

22 

27.  15 

39.  98 

42.  93 

63.89 

23 

7.55 

20.66 

7.49" 

26.60 

24 

9.  11 

28.46 

84.48 

9.  85 

25 

14.  38 

30.  29 

27.  97 

13.  01 

2 

X  critical  =  7,81«  Degrees  of  Freedom  =  3 
Indicates  non- rejection  of  differing  distributions 


» 


V 


21 


«  l  ”V  *%.  ■%.  ~w 


Table  2. 


C  Median  Test  Results  on  Season-to-Season  Differences 
n 


Altitude 

Barrow 

International  Falls 

Pt.  Mugu 

Brownsville 

1 

26.71 

67.  24 

16.  89 

47.  17 

2 

46.  73 

28.  65 

13.  05 

47.57 

3 

31.  15 

16.  66 

29.78 

28.55 

4 

8.27 

13.  75 

31.  29 

18.  08 

5 

21.64 

20.  22 

4.  03* 

11.  20 

6 

10.  21 

10.  53 

8.44 

12.64 

7 

39.  77 

6.  38 

9.  94 

18.  21 

8 

44.  38 

31.73 

7.  01 ' 

28.62 

9 

26.  68 

19.  45 

* 

5.69 

29.  31 

10 

9.82 

4.  97* 

11.  15 

15.  76 

11 

72.44 

8.  28 

59.39 

30.  37 

12 

81.95 

35.  29 

60.  03 

54.95 

13 

125.30 

92.43 

11.42 

30.  97 

14 

128.  77 

63.  91 

4.  74* 

26.59 

15 

188.  14 

74.44 

19.  93 

18.  22 

16 

207.  15 

118.09 

45.32 

28.  90 

17 

182.  60 

105.  92 

8.41 

78.  08 

18 

141.  02 

74.  05 

6.  61 

88.  15 

19 

144. 50 

65.  92 

3.63 

34.  73 

20 

188.  07 

94.  95 

23.  59 

4.49 

21 

161.49 

90.  86 

45.  84 

29.  70 

22 

136.62 

145. 68 

87.43 

70.  58 

23 

92.  28 

107.  82 

80.44 

51.  34 

24 

69.  79 

90.  39 

86.  52 

51.  62 

25 

107.43 

131.  05 

82.  37 

54.  10 

2 

x  critical  =  18#  Degrees  of  Freedom  =  3 


Indicates  non-  rejection  of  differing  distributions 


4.  THE  HIGH  ALTITUDE  DROP-OFF  RATE  FOR 


g 

Balsley  and  Peterson  first  observed  and  systematically  studied  the  fact  that 

o 

above  an  altitude  of  about  10  km  the  graph  of  the  logarithm  of  the  median  of  C 
versus  altitude  becomes  approximately  linear.  The  median  values  used  were 
obtained  from  radar  measurements  of  Cn  over  a  period  of  time  ranging  from  a  day 
to  several  weeks.  A  similar  observation  can  be  made  from  our  results  using 
Van  Zandt's  model  (see  Figures  1  to  8).  although  the  linear  region  appears  to  begin 
closer  to  15  km. 

2 

Defining  the  drop-off  rate  of  Cn  in  dB/km  by 

D-R  =  -  10/(z2  -  z x>  log  IC^(z2)/C^(ZjM  (18) 

where  z^  is  the  elevation  above  mean  sea  level  in  kilometers,  Balsley  and  Peterson 

then  look  at  the  variation  of  this  parameter  as  a  function  of  time  and  of  latitude. 

They  conclude  that  it  is  relatively  independent  of  time  and  is  a  decreasing  function 

of  latitude.  Based  on  their  somewhat  limited  data  they  suggest  that  much  more  work 

needs  to  be  done  to  clarify  these  conclusions.  Clearly,  if  valid,  they  go  a  long  way 

2 

towards  developing  a  climatology  for  Cn  in  the  upper  atmosphere  and  provide  a  way 
of  estimating  average  effects  of  the  turbulent  atmosphere  at  this  level  on  optical 
beams.  We  have  gone  on  to  investigate  this  area  extensively  by  use  of  Van  Zandt's 
model  of  1981.  5 

For  stations  of  various  latitudes  and  for  each  season  we  have  fitted  the  median 
2 

of  the  logarithm  of  CR  versus  z  to  a  straight  line  for  three  intervals:  (a)  10  to 
20  km  since  this  approximates  the  range  used  by  Balsley  and  Peterson;  (b)  15  to 
25  km  since  this  range  appears  to  be  the  best  straight  line  fit  to  the  Van  Zandt  model; 
and  (c)  over  the  entire  range  from  10  to  25  km.  Since  the  interval  from  15  to  25  km 
consistently  yielded  the  largest  value  of  r  ,  the  square  of  the  Pearson  correlation 
coefficient,  we  have  defined  our  drop-off  rate  for  this  interval.  See  Table  3  for 
typical  results  for  the  three  altitude  ranges  listed  above. 


2 

Table  3.  Drop-Off  Rate  (in  dB/km)  of  Cn  and  Associated  Square  of  Correlation 
Coefficient  for  Pt.  Mugu,  Calif. ,  1974,  for  Three  Altitude  Ranges 


Range 

W inte  r 

Spring 

Summer 

Fall 

10  -  20  km 

1.38  (0.993) 

1.  29  (0.  979) 

1.24  (0.965) 

1.35  (0.986) 

15-25  km 

1.52  (0.999) 

1.54  (0.997) 

1.54  (0.997) 

1.53  (0.999) 

10-25  km 

1.45  (0.  997> 

1.41  (0.991) 

1.37  (0.986) 

1.43  (0.994) 

We  first  calculated  the  drop-off  rate  for  six  stations  of  varying  latitude  for  the 
four  seasons.  Our  results  are  summarized  in  Table  4.  We  observe  little  seasonal 
variation,  especially  at  high  latitudes,  which  would  support  the  hypothesis  that  this 
parameter,  when  averaged  over  a  long  enough  time  period,  is  constant.  We  can 
also  see  that  it  is  a  decreasing  function  of  latitude  but  this  effect  is  not  as  pronounced 
as  it  is  in  the  radar  results  of  Balsley  and  Peterson.  This  could  be  due  to  the  fact 
that  our  averaging  is  over  a  much  longer  period  of  time,  about  90  days  for  each 
value.  In  Figure  9,  we  have  graphed  these  results  (as  well  as  two  additional 
stations,  at  about  48°  N  latitude,  see  below)  for  the  winter  seaso"  along  with  the 
radar  results.  We  chose  winter  because  most  of  the  radar  results  were  for  this 
season. 


9 

Table  4.  Drop-Off  Rate  (in  dB/km)  of  C  for  Four  Seasons  at  Various 
Latitudes,  197  4  n 


Station 

Latitude 

Winter 

Spring 

Summer 

Fall 

Barrow,  Alaska 

7  1°  28'  N 

1.  30 

1.  30 

1.  31 

1.  32 

International  Falls,  Minn. 

48°  34'  N 

1.  36 

1.40 

1.  44 

1.42 

Chatham,  Mass. 

41°  40'  N 

1.49 

1.54 

1.52 

1.48 

Pt.  Mugu,  Calif. 

34°  06'  N 

1.52 

1.  54 

1.54 

1.53 

Brownsville,  Tex. 

25°  54' N 

1.68 

1.69 

1.55 

1.63 

San  Andres,  Columbia 

12°  35'  N 

1.  71 

1.  58 

1.  59 

1.69 

Next,  we  calculated  the  season  drop-off  rates  for  three  stations  of  nearly  the 
same  latitude  but  differing  longitudes.  These  results  are  presented  in  Table  5. 
We  observe  little  difference  with  changing  longitude  and,  again,  little  seasonal 
variation. 
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Figure  9.  C„  vs  Latitude:  Radar  Results 
and  Van  Zandt  Model  Results  From 
Rawinsonde  Data  for  Winter  1974 


Table  5.  Drop-Off  Rate  (in  dB/km)  of  C  ,  1974,  for  Four  Seasons  at  Three 
Stations  of  Comparable  Latitude  n 


Station 

Latitude 

Longitude 

Winter 

Spring 

Summer 

Fall 

Sault  Sainte  Marie,  Mich. 

46°  28'  N 

84°  22'  W 

1.  37 

1.  43 

1.  46 

1.43 

Great  Falls,  Mont. 

47°  29'  N 

111°  21'W 

1.  37 

1.43 

1.48 

1.45 

International  Falls,  Minn. 

48°  34'  N 

93°  23'  W 

1.  36 

1.  40 

1.  44 

1.41 

We  can  also  determine  the  drop-off  rate  using  a  different  approach.  The 


models  of  Van  Zandt * 3  suggest  that  we  can  write 


c2  =  (P/T)2  G 
n 


where  the  function  G  depends  on  S,  N  and  L  (see  Section  3).  The  function  G  includes, 
among  other  things,  the  fraction  of  an  observed  layer  that  is  turbulent.  Since  P/T 
is  proportional  to  the  density,  which  we  know  decreases  with  increasing  altitude, 
it  is  of  interest  to  rewrite  the  above  as 


log  C2  =  log  (P/T)2  +  log  G 


and  calculate,  from  our  data  base,  the  density  drop-off  rate  (in  dB/km)  of  (P/T)  ; 
that  is. 


(D-R)  =  -10/<z2  -  Zj)  log[(P2/T2)2/(P1/T1)2] 


where  P.  and  T.  are  the  pressure  and  temperature  at  altitude  z..  This  parameter 
ii  2  1 

can  then  be  interpreted  as  that  part  of  the  drop-off  rate  in  due  only  to  decreasing 


density.  This  can  also  be  interpreted  as  the  drop-off  for  a  constant  CT-  For  the 


usual  exponential  atmosphere  with  a  constant  density  scale  height,  H^.  the  drop-off 
rate  is 


(D-R)  =  8.  69/  H  . 

P  P 


The  density  scale  height  is  related  to  the  tabulated  pressure  scale  height,  H^,  as 


H  =  H  /  1  +  H  In  T  . 


In  the  stratosphere  where  temperature  change  is  small,  11^=11^  -  6.4  km  and  the 
standard  atmosphere  drop-off  rate  is  (D  -  R>  =  1.  36  dB/km. 

P  2 

After  observing  from  our  data  base  that  log  (P/T)  is,  to  a  good  approximation. 


linear  from  z  =  15  km  to  z  =  25  km  we  evaluated  the  above  for  z^  =  15  km  and 


z2  =  25  km  for  the  same  six  stations  used  in  Table  4  for  each  season.  Since  very 
little  seasonal  variation  was  found,  in  Table  6,  we  present  just  the  average  values 


for  1974.  For  comparison  we  have  also  provided  the  average  value  of  the  drop-off 
2 

rate  in  C  for  1  " 4  in  Table  6.  This  table  demonstrates  that  the  two  drop-off  rates 


are  quite  close,  oiffering  by  less  than  1  percent  at  71°  N  and  by  about  7.  5  percent 


at  12°  N,  which  suggests  that  the  primary  cause  of  the  drop-off  rate  of  is  the  fact 


that  density  decreases  with  elevation.  We  further  note  a  decrease  in  the  drop-off 
2 

rate  of  (P/T)  with  increasing  latitude. 


I 


/.  *  »  »  , 
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Table  6.  Drop-Off  Rate  (in  dB/km)  of  (P/T)2  and  C2  for  Various 
Latitudes,  1974  n 


Station 

Latitude 

Drop-Off 

(P/T)2 

Drop^-Off 

cn 

Barrow,  Alaska 

71“  28'  N 

1.  30 

1.  31 

International  Falls,  Minn. 

48“  34'  N 

1.  37 

1.41 

Chatham,  Mass. 

41°  40'  N 

1.40 

1.51 

Pt.  Mugu,  Calif. 

34“  06' N 

1.44 

1.  53 

Brownsville,  Tex. 

25“  59'  N 

1.49 

1.  64 

San  Andres,  Columbia 

12“  35'  N 

1.  52 

1.  64 

5.  TRANSVERSE  COHERENCE  LENGTH 

The  transverse  coherence  length,  r  ,  is  a  measure  of  the  distortion  of  an 

optical  wavefront  by  atmospheric  turbulence.  Coherence  length  is  related  to  seven 

17 

specific  optical  properties  of  a  beam.  Over  a  circle  of  diameter  r  ,  wavefront 
distortion  has  an  rms  value  of  almost  exactly  1  radian.  The  coherence  length  is 
given  by  the  expression 

rQ  =  2.  1  (1.46  k2  I)*3/£  (24) 

where  k  is  the  wavenumber,  which  is  assumed  for  our  calculations  to  be  at  the  center 

2 

of  the  optical  band,  A  =  500  nm.  I  =  /  C  (z)Q  (z)  dz  where  Q  (z)  is  a  function  that 
depends  on  the  nature  of  the  optical  source.  For  an  infinite  plane  wave  source, 

Q(z)  =  1.  This  type  of  source  is  assumed  throughout;  then  in  order  to  calculate  rQ 
for  vertical  propagation  through  the  atmosphere  we  need  to  evaluate 

OO 

1  =  /  C'2  (z)  dz  (25) 

z  n 
o 

where  zQ  is  the  appropriate  altitude. 


17.  fried,  D.  L.  ,  and  Meyers,  G.  E.  (1974)  Evaluation  of  r  for  propagation  down 
through  the  atmosphere,  Appl.  Opt.  ,  13:2620-2622.  ° 


To  become  familiar  with  the  order  of  magnitude  of  the  quantities  involved,  we 

18 

refer  to  the  work  of  Gracheva  and  Gurvich  who  have  evaluated  this  integral  and 

then  rQ  under  the  extremes  of  strong  and  weak  turbulence.  For  these  cases  they 

obtained  an  I  of  1.  37  X  10  m*^  and  5.  79  X  10  ^  m*^  respectively.  Their 

respective  coherence  lengths  were  then  0.42  and  11.  1  cm  for  a  wavelength  of 

500  nm.  They  suggest  using  the  geometric  mean,  r  =  2.  2  cm,  as  a  representative 

1 9° 

value.  By  different  methods,  Walters  and  Kunkel  estimate  the  contribution  to  rQ 
of  this  integral  above  1  km  can  be  of  the  order  of  10  to  20  percent.  Thus,  for 
rQ  =  10  cm  at  2km,  the  contribution  above  1  km  is  between  1  and  2  cm.  Although 
these  values  may  be  helpful  in  order  of  magnitude  calculations,  in  many  applica¬ 
tions,  a  more  accurate  determination  of  rQ  is  required  to  reflect  local  conditions 
more  correctly. 

Although  we  expect  the  boundary  layer  to  make  the  major  contribution  to  this 

19 

integral  when  z^  lies  in  the  layer,  the  contribution  of  the  free  atmosphere  is  not 
negligible.  We  have  followed  two  approaches  in  calculating  this  contribution;  the 

5 

first  approach  is  based  on  Van  Zandt's  1981  model  and  the  second  is  based  on 

3 

Hufnagel’s  model. 

2 

By  use  of  the  Cn(z)  values  obtained  from  Van  Zandt's  model  (Section  3)  at  each 
km  level  from  z  =  1  km  to  z  =  25  km,  we  have  been  able  to  approximate  I  by  use  of 
Simpson1  s  Rule  over  the  range  from  zQ  =  1  km  to  zQ  =  20  km.  The  lower  value  of 
zQ  is  forced  on  us  due  to  the  lack  of  rawinsonde  data  below  1  km  (Section  3).  We 
stop  at  Zq  =  20  km  since  we  must  use  z  =  25  km,  not  infinity,  as  the  upper  limit  of 
integration  in  I.  One  can  estimate  that  this  approximation  is  accurate  to  about 
5  percent  at  z^  =  20  km  and  increases  rapidly  in  accuracy  for  zQ  <  20  km.  Should 
we  wish  to  extend  the  range  above  20  km,  we  could  easily  do  so  using  the  appropriate 
drop-off  rate  (Section  4).  This  is  probably  not  necessary  since,  at  this  high  an 
altitude,  Tq  is  so  large  that  diffraction  effects  dominate  over  optical  turbulence 
effects  in  determining  beam  distortion.  Some  typical  results  from  this  calculation 
are  shown  in  Figures  10  and  11. 


18.  Gracheva,  M. ,  and  Gurvich,  A.  (1980)  A  simple  model  for  calculation  of 

turbulence  noise  in  optical  systems.  Atmospheric  and  Oceanic  Physics, 
^6j819-822. 

19.  Walters,  D.  L. ,  and  Kunkel,  K.  E.  ( 198 1)  Atmospheric  transfer  function 

for  desert  and  mountain  locations:  the  atmospheric  effects  on  r  , 

J.  Opt.  Soc.  Am.,  Jl:397-405.  ° 
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Figure  10.  Coherence  L„ength  vs  Altitude:  Spring  1974,  Pt.  Mugu,  Calif 
Barrow,  Alaska;  Brownsville,  Tex.  ;  and  International  Falls,  Minn. 


COHERENCE  LENGTH  (CM) 


Figure  11.  Coherence  Length  vs  Altitude:  Winter  and  Summer  1974, 
Pt.  Mugu,  Calif. 


In  many  applications  an  accurate  determination  of  Tq  above  the  first  inversion 

layer  may  not  be  possible  due  to  lack  of  data  and/or  may  not  be  necessary  due  to 

the  limited  contribution  of  this  region  when  compared  to  that  of  the  region  below. 

19 

Walters  and  Kunkel  describe  an  approach  used  to  estimate  this  contribution.  We 
have  integrated  Hufnagel's  model  (Section  2)  to  evaluate  I,  obtaining: 


I  =  8.  2  X  10 


-23  2 


f(zn)  e 


+  4.  05  X  10' 


■13 


(-z0/1.5> 


(26) 


where 


This  should,  for  a  proper  choice  of  u,  form  a  reasonable  estimate  of  I  for  Zq 

corresponding  to  heights  above  the  first  inversion  layer. 

3 

Hufnagel  suggests  that  u  is  normally  distributed  and  for  Maryland  has  a  mean 
of  27  m/sec  with  a  standard  deviation  of  9  m/sec.  In  order  to  calculate  rQ  for 
conditions  of  high,  moderate,  and  low  turbulence,  we  have  used  the  above  expression 
for  I  with  u  =  36,  27,  and  18  m/sec  respectively.  The  results  are  shown  in 
Figure  12.  Although  the  graph  is  shown  down  to  zq  =  1  km,  it  must  be  remembered 
that  this  model's  validity  does  not  begin  until  we  are  above  the  first  inversion  layer 
which  may  or  may  not  be  this  low.  A  similar  calculation  can  be  done  at  other 
locations  using  the  appropriate  value  of  u. 


It  has  been  suggested  that  we  try  to  estimate  the  effects  of  the  atmosphere  below 

20 

the  first  inversion  layer  in  r  through  the  use  of  Kaimal's  model  as  modified  by 
19  °  21 

Walters  and  Kunkel  and  as  has  been  done  by  Brown  et  al.  This  unfortunately  does 
not  seem  feasible  using  our  data  base  since  the  rawinsonde  data  lacks  the  necessary 
resolution  required  to  consistently  determine  the  height  of  the  first  inversion  layer. 

In  order  to  compare  the  above  two  methods,  we  must  choose  a  method  for  deter¬ 
mining  a  seasonal  value  of  u  to  use  in  Hufnagel's  model.  At  each  kilometer  level, 
from  5  to  20  km,  for  each  launch  we  calculate  the  square  of  the  velocity.  We  then 
choose  the  median  for  that  level  for  the  appropriate  season.  It  is  felt  that  because 
of  the  spread  in  values  of  wind  speeds  from  launch  to  launch,  the  median  is  more 
representative  than  the  mean.  Then  u  is  the  rms  value  of  these  medians  from 
5  to  20  km. 

2 

Agreement  between  the  two  models  is  not  particularly  good.  When  u  lies  in 
2  2 

the  range  150  to  275  m  /sec  agreement  is  best  (Figure  13).  At  higher  values  of 
2 

u  ,  Van  Zandt's  model  yields  a  value  of  r^  in  excess  of  Hufnagel's  model  for  the 

same  value  of  z^  (Figure  14)  while  at  low  wind  speeds  the  reverse  is  true  (Figure  15). 

The  reason  for  this  discrepancy  is  that  Van  Zandt's  model  is  relatively  insensitive 
2  2 

to  u  whereas  for  values  of  z  and  u  sufficiently  large,  Hufnagel's  model  predicts 
2 

that  Cn  should  be  directly  proportional  to  this  parameter.  Figure  16  is  a  graph  of 
versus  z  for  Chatham,  Mass,  for  the  winter  of  1974  and  San  Andres,  Columbia 
for  the  spring  of  1974.  These  two  stations  and  seasons  were  chosen  for  comparison 

p 

because  they  have  the  largest  and  smallest  values  of  u  that  we  have  encountered. 

2  2  2  2 

1258  m  /sec  and  36.6  m  / sec  respectively.  We  see  little  difference  between  the 
two,  but  Hufnagel's  model  predicts  that  they  should  differ  by  a  factor  of  34. 


6.  CONCLUSIONS 


We  have  used  Van  Zandt's  model  of  1981  (with  some  modifications)  in  con- 

junction  with  the  rawinsonde  data  base  to  calculate  median  values  of  C  at  various 

stations  for  each  season  at  each  kilometer  level  from  1  km  to  25  km.  These 

results  were  examined  statistically  for  seasonal  and  latitudinal  variations.  We  have 

2 

calculated  the  drop-off  rate  of  C  (in  dB/km)  from  15  km  to  25  km  and  the  coherence 

n 

length,  rn,  for  plane  wave  propagation  down  through  the  atmosphere.  In  addition, 
u  13 

we  used  Hufnagel’s  model  to  calculate  r^  from  the  same  database.  We  conclude 
that: 


V  1 


20.  Kaimat,  J.C.  Wyngaard,  J.C. .  Haughen,  D.  A.,  Cote,  O.  R.  ,  Izumi,  Y.  L.  , 

Caughey,  S.  J. ,  and  Readings,  C.J,  (197  6)  Turbulence  structure  in  the 

convective  boundary  layer,  J.  Atm.  Sci. ,  33:2152-2169. 

-  v*. 

21.  Good,  R.  E.  ,  Brown,  J,  H. ,  and  Quesada,  A.  F.  ( 1982)  Measurement  of  high 

altitude  resolution  profiles  and  their  importance  on  coherence  lengths. 
Proceedings  of  SPIE,  the  Inter.  Soc.  for  Opt.  Eng.  365:105-111. 
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COHERENCE  LENGTH  (CM) 


Figure  13.  Coherence  Length  vs  Altitude:  Spring  1974,  Pt.  Mugu,  Calif. 
From  Van  Zandt  and  Hufnagel  Models 


2 

(1)  There  are  statistically  significant  differences  in  median  C^ 
distributions,  both  latitudinally  and  seasonally.  The 
four  stations  were  broken  up  seasonally  and  the  four  stations 
were  compared  for  each  season.  Each  station  was  also 
examined  separately  and  the  four  seasons  were  compared 
for  each  station. 


U**2  (MED I AN-OF-SQUARED- WI NDS  ) 


707.39 


COHERENCE  LENGTH  (CM) 


Figure  14.  Coherence  Length  vs  Altitude:  Winter  1974, 
international  Falls,  Minn. ,  From  Van  Zandt  and  Hufnagel 
Models 

(2)  The  drop-off  rate  is  a  decreasing  function  of  latitude,  but 
varies  only  slightly  with  longitude  and  season.  Seasonal 
values  from  Van  Zandt' s  model  are  somewhat  smaller  than 
those  obtained  from  radar  at  comparable  latitudes.  This  is 
possibly  due  to  the  shorter  time  periods  used  in  the  radar  results. 

(3)  Much  of  the  drop-off  rate  can  be  explained  as  due  to  the  decrease  in 
density  with  altitude,  99  percent  at  71*  N  and  about  92.  5  percent  at  12°  N. 


ALTITUDE  (KM) 


U**2  (MED 1  AN -OF -SQUARED- WINDS  )  =  55.77 


Figure  15.  Coherence  Length  vs  Altitude;  Summer  1974, 
Brownsville,  Tex.  ,  From  Van  Zandt  and  Hufnagel  Models 


(4)  The  coherence  length  for  an  optical  plane  wave  can  be  calculated 
from  Van  Zandt’s  model  and  from  Hufnagel's  model.  However, 
major  discrepancies  between  the  two  models  occur  at  very  low 
and  high  wind  speeds. 
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SUBROUTINE  VZCNSQ(NBAR,SBAR,CN2) 

DIMENSION  L(1?),N(100f ),S(1ff0) 

DIMENSION  PSt tfff) ,PN< 1000 > 

REAL  L,NBAR,N,MMBSIt 

DATA  L/0.2,0.4,0.4,0*8,1.0,2.t,3.t,4.f,A.f,10.,20.,40.,40.,80., 
XI  f f. ,2ff .  ,400. , 60 f . , 800 ./ 

FUNCTIONS 

C3  REPRESENTS  THE  THICKNESS,  L 

A  REPRESENTS  THE  STABILITY,  N 

ALL  DIMENSIONS  MUST  BE  IN  H,K,S,  SYSTEM 

SBAR  REPRESENTS  THE  AVERAGE  SHEAR  SQUARED 

H=-A  P/T  N/6 

6.266E-11  =  <A/G)**2 

CN  IS  REALLY  CN  <T/P) 

CN  IS  FOR  DRY  CONDITONS 


—  HERE  DEFINE  LOCAL  FUNCTIONS  SIGMAS  AND  CNSQ 
SIGMAS(ABAR,C3)=0.2*ABARm0.25*(C3>**<-0.3> 
CNSO<A,C3)=2.8*C3**<1.3333)*<6.266£-11>*A*A 

SQRT2PI=SQRT<2.0»3. 14159) 

RI=f .25 
SUM=f.f 
NCNT=20 
TSTART=-3. 

SQRTSBA=SQRT(SBAR) 

DO  4f  KL=2, 1 4 

THIS  LOOF-  INTEGRATES  OVER  THE  VARIOUS  SPECIFIED  LAYER  THICKNESS 
PL=0.01 

IFtKL.EQ.I)  THEN 
DEL-0 .2 
GO  TO  6 
END  IF 

DEL=L(KL)-L(KL-1 ) 

6  SIG=SIGMAS(NBAR,L(KL) ) 

SI6SQ-SI6*SIG 
SIGMAN=SIG*SQRI (NBAR) 

SUHN=f .0 
DO  30  1=1 ,NCNT 

THIS  LOOP  INTEGRATES  OVER  THE  STABILITY  PARAMETER 


IFd.NE.I )  GO  TO  14 


DELT  =  -2  .•♦TSTAR  T/NCNT 

DELN=D£LT*SIGHAN 

T=TSTART 

N(I)=T*SI6f1AN+NBAR 
60  TO  IS 
T=T+BELT 

IF(T.GT.-TSTART)  GO  TO  31 
N( I )=T*SI6MAN+NBAR 
C2-T*T/2 

IF(N<I1.LT.-.0004)  GO  TO  30 
C»-C«SQ(N(I!,L(KL)) 

IF (C2.GT.1 10 . )  THEN 

PN(I)=0.0 

GO  TO  30 

ELSE 

PN(I)=EXP(-C2)/$QRT2PI/$I6MAN 
END  IF 

SUHS=1.0E-30 
DO  20  J~1 ,NCNT 

THIS  LOOP  IS  OVER  THE  S  VARIABLE 

IF(J.NE.I)  GO  TO  V 
S( 1 ) =N< I 1/R1 

IF(N(I) .LE.0.0)  S(1 )=1 .IE-8 
TT*(SQRT(S<1 ))-SQRTSBA)/SIG 
DELTT=-2.0*TSTART/NCNT 
C2=2.0»$IG»DELTT 
R=M<1 >/S<  J) 

GO  TO  18 
TT*TT+DELTT 

S<J)  =  <SI6*mSQRTSBA>+*2 
R=N( I )/S( J) 

IFfTT.GT.-TSTART. AND. R.LT. 0.05)  GO  TO  21 

IF (R. GT ,25 )  GO  TO  20 

DELS- ADS ( C2» ( SIG+TT*SQRTSBA) ) 

SS=SQRT(S( J)»SBAR)/SIGSQ 

iF( SS.GT. 473. )  GO  TO  20 

Cl =( SORT ( S( J) ) -SQRTSBA) «*2/(2.0*SIGSQ) 

IF<C1 -GT.100. )  THEN 

PS(J)=0.0 

GO  TO  20 

ELSE 

BESSEL=NHBSI0(2,SS,IER) 

THE  ARGUHENT  2  REPRESENTS  EXP(-SS)  KSS) 
PS(J)=EXP(-C1 )*BESSEL/(2.0*SIGSO) 

END  IF 

SUHSS=PS( J)*PL*CN*PN( I ) ♦DELS*DELN*HEL 
SUNS- SUNS +SUNSS 

IFINUI.GE. 0.0. AND. SUHSS/SUMS.LT. 0.001)  GO  TO  21 

CONTINUE 

SUHN-SUNN+SUHS 


C  IFISUHS'SUNN. IT. 0.001 >  GO  TO  31 

30  CONTINUE 

31  SUH=SUN+SUHN 

IF (SUNN/SUN. LT. 0.101)  GO  TO  41 

40  CONTINUE 

41  CN2=SUH 
RETURN 

999  FOFtHATU0X,3HKL=,14,8613.3) 

998  FORHAT (8X,3H  I=, 14,8015.3) 

99?  F0RNAT(2X,3H  J=,I4,8G14.3> 

996  FORHAT (2IX, 'SBAR  AND  NBAR  =  ' , 61 5. 4, 2X ,G1 5.6 ) 

995  FCRHAT120X,'  CN  SQ  =',015. 4) 

994  F0RHAT(A6) 

993  FORHAT (1  HI) 

992  FORHAT ( 1 IX,20H  TSTART  AND  NCNT  =  ,F6.0,I5) 

END 


COHHAND- 


